Electronic Inhomogeneity and Vortex Disorder in Superconducting Sr .75K .25Fe2As2 
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We use scanning tunneling microscopy to investigate the surface structure, superconducting and 
vortex properties in the hole-doped superconductor Sro.75Ko.25Fe2As2 (T c =32 K). We resolve the 
long debate over the nature of the ^4Fe2As2 cleaved surface, finding a dominant Sr/K termination 
with 1x2 reconstruction and ubiquitous superconducting gap, with rarer patches of gapless, 
unreconstructed As termination. The superconducting gap varies by a/ A= 16% on a ~3 nm length 
scale, with average 2A//c j bT c =3.6 in the weak coupling limit. At 9 T we observe isotropic vortices 
which provide a measure of the superconducting coherence length £ = 2.8 nm, and constrain the 
pairing symmetry. Furthermore, we observe a vitreous vortex phase and quantify its correlation 
length in comparison to other iron-based superconductors. The comparison leads us to suggest the 
importance of dopant size mismatch as a cause of dopant clustering and a new avenue to optimize 
vortex pinning for increased critical current. 

PACS numbers: 68.37.Ef, 74.55,+v, 74.70.Xa, 74.25.Uv 



The recent discovery of high transition temperature 
(T c ) iron-based superconductors (Fe-SCs) [1] has 
provoked much excitement in condensed matter physics, 
and launched a new era in the search for the key to 
high-T c superconductivity Q. Like cuprates, Fe-SCs 
exhibit a layered structure with electronically active 
superconducting planes separated by buffer layers, and 
the superconductivity develops from antiferromagnetic 
parent compounds upon chemical doping. In addition 
to enabling superconductivity, the dopants are potential 
sources of nanoscale phase separation Q , and crystalline 
[1] and electronic disorder [EH?}; which may in turn lead 
to T c suppression [8| and vortex pinning [9[. With the 
diversity of possible dopants in Fe-SCs, it has remained 
an elusive challenge to characterize and categorize their 
nanoscale effects. 

The scanning tunneling microscope (STM) is an ideal 
tool to study the nanoscale properties of correlated 
electron materials. STM studies of Fe-SCs have 
presented several controversial results [To[. First, the 
cleaved AFe2As2(001) surface showed both 1x2 and 
y/2 x y/2 reconstructions with unclear origin: either a 
1/2-layer of A or a reconstruction of a complete As layer. 
Second, spectroscopic images of optimally electron-doped 
Ba(Fei_ a; Co ;z; )2As2 cleaved at ~ 25 K revealed nanoscale 
variations in the superconducting gap A on a length 
scale of several nanometers [5j. However, studies of 
the same compound cleaved at room temperature found 
a shorter A correlation length of ~0.8-1.0 nm, closely 
matching the average Co separation for a random dopant 
distribution (~1.0 nm), which prompted the hypothesis 
that gap variations may be caused by the disorder of 
individual Co atoms [7]. Finally, the 1x2 surface of 



electron-doped BaFe1.8Coo.2As2 displayed a disordered 
vortex lattice without observed vortex core Andreev 
bound states Q. In contrast, the domain-like surface 
of hole-doped Bao.6Ko.4Fe2As2 displayed a hexagonal 
vortex lattice with pronounced vortex core bound states 



U| . It has been proposed but not verified that the vortex 



discrepancy may be explained by stronger scattering from 
the in-plane Co dopants than the out-of-plane K dopants. 

To address these controversies, Sri_ :c K iC Fe2As2 is a 
unique system with specific advantages. Gao et al 
predicted no surface reconstruction in the As-terminated 
SrFe2As2 surface [12[ (in contrast to As-terminated 
BaFe2As2). This makes it easy to distinguish between 
a complete As layer and a partial Sr layer. Moreover, 
gap and vortex mapping in this material can test 
the hypothesis that individual strong-scattering Co 
dopants are responsible for the gap inhomogeneity and 
vortex pinning differences between BaFe1.8Coo.2As2 and 
Ba .6Ko.4Fe 2 As2. 

Here we utilize a home-built cryogenic STM to 
systematically study the hole-doped superconductor 
Sro.75Ko.25Fe2As2 (underdoped, T c =32 K). The single 
crystal samples we re p repared by a FeAs flux method 
detailed elsewhere [13|- We cleave the samples in situ 
at ~25 K and insert them immediately into the STM 
for imaging at 6 K. To obtain a tunneling current, we 
apply a bias to the sample while the tip is held at 
virtual ground. All spectra were measured using a lock-in 
technique with a 1.0 mV rms bias modulation at 1110 
Hz. Figure 1(a) shows the schematic crystal structure 
of Sri_ ;E K ;E Fe2As2. Due to the strong covalent bond 
between Fe and As ions, the cleavage likely occurs in 
the Sr/K plane, leaving ~l/2 Sr/K on either exposed 
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side to balance the chemical valence. Figure 1(b) typifies 
the most commonly observed surface with local 1x2 
stripes [3, [H|. Occasionally for low-T cleaves, the 
metastable As-terminated lxl surface may be expected 
[l2| . Previous STM imaging of SrFe2As2 samples cleaved 
at 77 K showed no such lxl patches [14], and instead 
both 1x2 and a/2 x y/2 orderings in the entire cleaved 
surface were explained as the bare but reconstructed 
As layer; whereas a previous study of Sri_ x K :c Fe2As2 
cleaved at 10 K showed only a very small lxl patch 
with no accompanying spectroscopy [15j. However, we 
observe larger lxl patches [Fig. 1(c)], constituting 
~ 5% of the surface, with average As- As atomic spacing 
of approximately 0.4 nm. This first observation of large 
areas of 1 x 1 ordering provides the final unambiguous 
proof that the dominant 1x2 and y/2 x y/2 observed 
structures arise from a 1/2-Sr/K layer. The bright rows 
in our lxl regions are identified as residual Sr/K atoms. 
We then record differential dl/dV spectra on both Sr/K- 
and As-terminated surfaces, illustrated in Figs. 1(d) and 
1(e), respectively. In stark contrast to the universal 
superconducting gap with clear coherence peaks on the 
Sr/K-terminated surface, no symmetric superconducting 
gap is observed on the As-terminated surface, which 
is probably caused by the strong polarity of the latter 
surface. The cleaved structure therefore plays a crucial 
role in superconductivity at the surface, in some cases 
preventing even proximity-induced superconductivity 
from appearing. 

To further explore the effect of cleaved surface 
structure on superconductivity, we survey dl/dV spectra 
at all locations within the cyan box in Fig. 1(b). 
These spectra can be analyzed to yield maps of A(r), 
zero bias conductance (ZBC) Z(r), and coherence peak 
strength C(r), as displayed in Figs. 2(a)-2(c). All 
three maps show spatial inhomogeneity. The average 
A = 5.0 meV with a standard deviation of <r=0.8 meV 
gives a fractional variation a/A = 16%, larger than the 
12% variation found in BaFe1.8Coo.2As2 [5[. The small 
reduced gap 2A//cbT c ~3.6 suggests weak coupling in 
Sro.75Ko.25Fe 2 As 2 . 

In order to investigate the spectral shape, all spectra 
are binned based on A, with their averages plotted in Fig. 
2(d). Evidently, the spectra with smaller A tend to show 
larger ZBC and stronger coherence peaks. To quantify 
this trend, we compute the cross-correlation between 
A(f) and Z(r) as well as coherence peak strength C(r) 
[Fig. 2(e)]. Together with the autocorrelation of the 
gap map, we note that the characteristic length scale 
for all correlations is ~3 nm, exceeding the average 
separation (~0.8 nm) between individual K dopants. We 
therefore suggest that the K dopants exhibit nanoscale 
phase separation fl6l |. clustering to form K-rich and 
K-poor regions. This can more effectively relax the 
strain caused by the large ion size mismatch between 
K+ (146 pm) and Sr 2+ (126 pm) 0. Due to the 
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FIG. 1. (color online) (a) Schematic crystal structure of 
Sri_ x K x Fe 2 As 2 . (b) STM topography (V s = -100 mV, I = 
35 pA) of the commonly observed Sro.75Ko.25Fe2As2 surface 
with 1x2 stripe order. Inset shows a zoom-in of the stripe 
with the magenta dots denoting Sr/K atoms ( V s = -100 mV, I 
= 30 pA, 2 nm x 2 nm). (c) STM topography ( V s = -100 mV, 
/ = 350 pA) showing As-terminated lxl surface, decorated 
by sparse Sr/K rows. Inset shows magnification of the lxl 
surface with the yellow dots denoting As atoms (V s = -100 
mV, I = 350 pA, 2 nm x 2 nm). (d, e) Spatially averaged 
dl/dV spectra in (b) and (c) regions, respectively. Tunneling 
gap was stabilized at V s — -100 mV and / = 300 pA. 



high sensitivity of superconductivity on the chemical 
doping [l8|, the resultant K clusters could account well 
for our observed electronic inhomogeneity. Meanwhile, 
the apparent anticorrelation between A(r) and C(f) may 
stem from the presence of the competing SDW gap 
[Hi EM I n regions of smaller K-doping, the local 

SDW gap would be relatively enhanced, which would 
suppress Z(r) and reduce the low energy states available 
to contribute to C(r), and furthermore cause these 
reduced superconducting coherence peaks to appear as 
shoulders shifted outwards on the larger- SDW-gapped 
background, enhancing A(r). 

Magnetic vortices are technologically important, as 
the superconducting critical current J c is limited by the 
vortex pinning strength. Vortices are also scientifically 
valuable for determining the superconducting coherence 
length £ 0, [ll| and pairing symmetr y lioj , and serving 
as windows into competing phases [21[. Figure 3(a) 
shows the vortices in a 9 T c-axis magnetic field as 
a map of dl/dV at the filled state coherence peak (-5 
meV). The vortices locally suppress the superconducting 
coherence peaks, and appear as purple-black features 
with depressed conductance. To better emphasize the 
vortices, Voronoi cells are overlaid onto the image. From 
the Voronoi cell size, we estimate the average flux per 
vortex = 2.1 ± 0.1 x 10 -15 Wb, which is consistent 
with one magnetic flux quantum, <I>o = 2.07 x 10 -15 Wb. 

In sharp contrast to the hexagonal vortex lattice 
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FIG. 2. (color online) (a-c) Maps of A(f) (half the 
distance between coherence peaks), Z(f), and C(f) (average 
conductance at the two peaks), (d) Binned and averaged 
spectra for five ranges of gap A (A m i n -3.8, 3.8-4.7, 4.7-5.5, 
5.5-6.3, 6.3-A m ax meV from top to bottom), color-coded to 
match those in (a). The spectra have been vertically offset 
for clarity, with their zero-conductance positions marked by 
correspondingly colored horizontal lines. Vertical dashes at 
±5 meV are guides to the eye. (e) Azimuthally averaged 
autocorrelation of A(r), and cross-correlations of A(r) with 
Z(f) and C(f). Raw spectra are used in (d), but spectra were 
normalized by their backgrounds to compute Z(f) and C(f) 
in (b) and (c). Backgrounds were extracted from a cubic 
fit to the conductance for |V| > lOmV. The normalization 
procedure does not affect the signs of the correlations in (e). 



which indicates negligible pinning in Bao.6Ko.4Fe2As2 
ll|, the vortices in Sro.75Ko.25Fe2As2 do not form 
an ordered lattice. Instead, similar to electron-doped 
Ba(Fei_ :c Co :c )2As2 [H, [22|, short-range hexagonal order 
(vortex glass phase) is justified based on the following 
two tests. First, the relative distances da = |r 7 - — 



rj\ are calculated for all vortex pairs at positions 
Ti and rj. Figure 3(b) plots the histogram of the 
observed dij. The pronounced peak at the smallest 
distance of ~16.4 nm coincides almost exactly with 



the expected lattice constant a = y 2<&o/ \/3H (~16.3 
nm, the first vertical line) for a perfect hexagonal 
vortex arrangement at 9 T. However, with increasing 
distance the experimental histogram peaks cease to 
match the corresponding hexagonal lattice separations. 
This presents the first indication that the vortices are of 
short-range hexagonal order. Second, Fig. 3(c) shows 
the Delaunay triangulation by connecting all nearest 
neighbor vortex sites. Each vortex, surrounded by a 
closed Voronoi polygon, is color coded in terms of the 
coordination number. The statistics of the coordination 
number and angles of Delaunay triangles are plotted in 
Figs. 3(d) and 3(e), respectively. Two-thirds of vortices 
are six-fold coordinated, and more importantly, the 
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FIG. 3. (color online) A 80 nm x 80 nm dl/dV map at V s = 
-5 meV, recorded in a 9 T magnetic field. The vortex centers 
are determined by 2D Gaussian fitting, (b) Histograms of 
the vortex pair distances dij from (a) and another similar 100 
nm x 100 nm vortex map. Vertical blue lines correspond to 
the positions and relative weight of pair distances for an ideal 
hexagonal vortex lattice at 9 T. (c) Delaunay triangulation 
(black lines) of the vortex lattice in (a). Each vortex is 
color-coded based on the number of its nearest neighbors, (d, 
e) Coordination number and Delaunay angle distributions on 
the border- free regions. Solid blue lines are the Gaussian fits. 



Delaunay angle distribution shows a single pronounced 
peak at ~ 60°. All this evidence consistently supports a 
short-range hexagonal order of the vortex arrangement, 
indicative of strong vortex pinning in Sro.75Ko.25Fe2As2- 

To gain insight into the internal vortex structure, 
we register all vortex centers, then average the density 
of states around 48 vortices, as depicted in Fig. 
4(a). In contrast to the two- and four- fold symmetric 
vortices in FeSe [20] and LiFeAs [23j, the vortices in 
Sro.75Ko.25Fe2As2 are nearly isotropic. This observation 
does not support the claims of d- wave pairing in more 
overdoped 74i_ :E K ;E Fe2As2 materials [24] . although it 
leaves open the possibility that the isotropy stems 
from thermal smearing [23| or impurity scattering (25| . 
Figure 4(b) shows a series of differential conductance 
dl/dV spectra across one vortex. As expected, Z(r) is 
elevated within the vortex core, reflecting the suppression 
of superconductivity there. Figure 4(c) plots the 
vortex-induced Z(f) as a function of the radial distance 
from the vortex centers. An exponential decay of Z(r) = 
Zoo + Aexp(— r/£) fits the data well, resulting in a 
superconducting coherence length of£ = 2.8±0.5 nm. 

Now we shed some light on the vortex pinning in 
Sro.75Ko.25Fe2As2. Charge doping has recently been 
argued to lead to collective pinning of vortices in 
Fe-SCs (26[. However, this charge model alone cannot 
account for the contrast between our observed vortex 
pinning in Sro.75Ko.25Fe2As2 and the ordered vortices 
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TABLE I. Ion size and vortex arrangement for several Fe-SCs. 
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FIG. 4. (color online) (a) Average dl/dV at -5 meV from 
overlaying 48 single vortices, (b) Typical series of normalized 
dl/dV spectra straddling a single vortex. The spectra are 
equally separated and span a total distance of 11 nm. (c) 
Radial dependence of Z(r) around three vortices. 



in Bao.6Ko.4Fe2As2 with identical K dopants. Another 
possibility is that the shear modulus Cqq, rou ghly 
proportional to H c2 b(l - b) 2 (where b = H/Hno)\21\. 
may play a role in the different vortex arrangements 
between Bao.6Ko.4Fe2As2 and Sro.75Ko.25Fe2As2. A 
larger Cqq often results in an ordered vortex lattice. 
Using the Ginzburg-Landau expression H C 2 = <3>o/27r£ 2 , 
we estimate H C 2 = 42 ± 15 T in Sro.75Ko.25Fe2As2, 
whereas H c2 ~ 75 T in Ba .6Ko.4Fe 2 As2 [n|. The 
persistence of vortex lattice order in Bao.6Ko.4Fe2As2 
down to 4 T (where the smaller Cqq is comparable 
to that of Sro.75Ko.25Fe2As2 at 9 T) [11] suggests 
that Cqq can not solely account for the short-range 
vortex order in Sr .75K .25Fe2As2. Finally, we consider 
the nanoscale electronic inhomogeneities caused by K 
clustering, as hypothesized earlier. Note that the 
gap correlation length of ~3 nm is comparable to the 
size of a single vortex core (£ ~ 2.8 nm). This 
indicates that the observed electronic irregularities may 
pin vortices [28| and disorder the vortex lattice. Since 
the ion size mismatch between K + and Ba 2+ is five 
times smaller than that between K + and Sr 2+ (see 
Table [I]), one can expect K + ions to be less clustered 
in Bao.6Ko.4Fe2As2 than in Sro.75Ko.25Fe2As2 [29|. This 
leads to greater homogeneity and weaker vortex pinning, 
consistent with the well-ordered vortex lattices in 
Bao.6Ko.4Fe2As2 [11]. In contrast, the large ion size 
mismatch between K + and Sr 2+ could account for the 
greater electronic inhomogeneity and disordered vortex 
lattice in Sro.75Ko.25Fe2As2 observed here. 

To quantify the effect of the ion size mismatch, we 
compute the correlation length of the vortex lattice ( 
=39 nm in Fig. 3(b) and several other Fe-SCs for 
comparison in Table [I] [22[. The larger ( = 128 nm in 
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Bao.6Ko.4Fe2As2 [Hi can be explained by the smaller 
ion size mismatch, whereas the small ( = 22 nm in 
BaFe1.8Coo.2As2 [5] can be explained by the large ion 
size mismatch there. Compared to Sro.75Ko.25Fe2As2 
which has a slightly larger ion size mismatch, the smaller 
( in BaFe1.8Coo.2As2 suggests that Co 3+ themselves, 
doped directly in the superconducting FeAs layer, may 
act as stronger pinning sites. Our hypothesis about 
the importance of dopant size mismatch is further 
supported in the FeSe^Tei-^ system, where there is a 
large size mismatch between Se 2_ and Te 2_ . Although 
a homogeneous superconducting gap and ordered vortex 
lattice have been demonstrated in stoichiometric FeSe 
[2o| . both nanoscale chemical phase separation (30| 
and a disordered vortex arrangement [3l| are seen in 
FeSe^Tei^. 

Finally we comment on the absence of quasiparticle 
bound states within the vortex cores [Fig. 4(b)] (32| . 
which often appear as a pronounced peak at or near Ep 
in other Fe-SCs 0, Ej . Using the residual resistivity 
po = 0.08 m£l • cm and Hall coefficient Ryl = 1.16 x 10" 9 
m 3 /C [HI, 0], we obtain the electronic mean free path 
£ = fi(37r 2 ) 1 / 3 /e 2 n 2 / 3 po °° 5.2 nm, two times bigger 
than £ ~ 2.8 nm. This suggests that our sample is 
macroscopically in the clean limit, where the vortex 
core bound states should have been observed. However, 
the vortices may be pinned in the relatively disordered 
regions, where the local mean free path is smaller (26| . 
Moreove, we note that the visibility of vortex core bound 
states apparently depends on the surface structure within 
the same material Bao.6Ko.4Fe2As2 [Til l33|. 

In summary, our detailed STM/STS study of surface 
structure, superconducting gap, and vortex arrangement 
in Sro.75Ko.25Fe2As2 has addressed four important 
questions on Fe-SCs. First, we have shown images and 
spectroscopy of the first large patches of unreconstructed 
lxl surface on any AFe2As2 material. This observation 
provides the final unambiguous evidence that both 1 x 
2 and x a/2 reconstructions seen before represent 
a partial Sr layer. We observe no superconducting 
gap on the 1 x 1 As surfaces, in contrast to the 
ubiquitous gap on 1 x 2 surfaces, which reiterates the 
importance of attention to surface details when using 
STM/STS to study bulk superconductors. Second, we 
have explicitly shown by spatially resolved spectroscopy 
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that the gap varies on a 3 nm length scale, quite 
different from the 0.8 nm average distance between 
individual K dopants. This supports a local K clustering 
model. Third, we have imaged a vitreous vortex 
phase with a short-range hexagonal order (£ ~ 39 
nm), suggesting intermediate pinning strength between 
BaFe1.8Coo.2As2 and Bao.6Ko.4Fe2As2- Our hypothesis 
of the importance of dopant size mismatch provides a 
new avenue to optimize vortex pinning for increased 
critical current J c . Fourth, we demonstrate an isotropic 
vortex structure and superconducting coherence length, 
£ ~ 2.8 nm. This constrains the pairing symmetry in 
Sro.75Ko.25Fe2As2, and helps to understand the nature 
of the superconducting mechanism in Fe-SCs. 
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